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FOREWORD 


The material documented in this interim report describes 
progress to date accomplished under NASA Contract No. NASS- 
319S0f during the period 13 April 1976 through 14 June 1977. 
During this period, a methodology has L^en formulated, and 
a general computer code implemented and checked out, for 
processing sinusoidal vibration test data to simultaneously 
make adjustments to a prior mathematical model of a large 
structural system, and resolve measured response data to 
obtain a set of orthogonal modes representative of the test 
article. 

The general procedure is referred to as " model verification ". 
The term "model verification" is used herein to denote a .. 
procedure with two distinct and equally important objectives: 
(1) to establish a proper model configuration by examining 
different variations of configuration with respect to their 
ability to match available test data, and (2) to estimate 
specific parameter values for a selected model configuration. 
The first objective is met by providing a general modeling 
capability within the logical structure of the computer code. 
The practical utility provided by this modeling capability 
is intended to facilitate a "man-in-the-loop" type of function 
where the user may apply his judgement and modeling skill to 
achieve a proper model configuration. The second objective 
is met by providing fully automated parameter estimation 
programs to optimize the fit of any selected model configura- 
tion to the given test data. 

Consistent with the "nuui-in-the-loop" approach, the computer 
code is implemented in two separate programs, where the out- 
put from one becomes input to the other. Thus, the user may 
inspect output from the first (frequencies, damping, mode 


shapes and sinusoidal response) prior to completing computa- 
tions to evaluate revised mass and stiffness matrices# and 
the changes which have been made to the prior mathematical 
(matrix) model. 

During the continuing part of this contract# modeling and 
computational procedures will be refined to improve the 
overall capability of the procedure to deal with real«world 
problems# and to provide user-oriented guidelines for 
application of the methodology. 
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INTRODUCTION 


The J. H. Wiggins Company has developed, under contract to 
NASA-MSFC, a methodology auid computer program £or structural 
system identification. This work addresses a general need 
in the area of analysis-test correlation with particular focus 
on the Space Shuttle Program. The objective of this effort is 
to enable computerized processing of sinusoidal vibration test 
data to refine pre-established analytical models. The two 
primary objectives of the effort ares (1) to develop and check 
out the computer code itself; and (2) to demonstrate its 
application to a real problem, i.e., using actual test data 
cmd a mathematical model for a real structural system* 

The first objective has been achieved. A computer code has 
been developed emd checked out using simple analytical models 
and simulated test data* Several analytical models have been 
investigated in an effort to gain basic experience with the 
parameter estimation procedure, and to examine numerical 
sensitivity to different types of parameters. This experience 
has verified the basic operation of the code. In addition, 
it has pointed to some areas where improvement to the existing 
code is desirable from the standpoint of practical applica- 
tion to real structural systems. 

The proposed demonstration problem has not been undertaken. 

It was originally proposed that the full-scale orbiter and 
HGVT data be used. Although the present computer code can 
be applied to this problem, it became apparent that the over- 
all complexity of the problem, considering the dynamic model 
as well as the HGVT data, is such that the likelihood of 
achieving the second objective is not as good as originally 
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hoped. Consequently a different tack is proposed: Select a 

simpler demonstration problem first and gradually increase 
problem complexity until the full-scale orbiter data can be 
evaluated. 

Three basic steps are involved in the present computerized 
methodology with an additional step being recommended here. 

1) Establish a suitable prior models including both a 
modal representation (the analytic model) and mass 
and stiffness representation in terms of the physi- 
cal coordinates (the dynamic model) . 

(NEW) perform a "zero" order identification using measured 
natural frequency and mode shape data. 

2} Refine the parameters of the analytic model using 
measured response data. 

3) Refine the parameters of the dynamic model using the 
refined amalytic model as if it were test data. 

It has been determined that the "closeness” of the prior model 
to the actual system is an important factor in the convergence 
behavior of the estimator. This is especially true in the case 
of the Orbiter where response data are only available at 
resonant frequencies. Such data will constitute a rather 
sparce sampling of frequency response characteristics, tending 
to make convergence more difficult to achieve. 

It has been determined that some additional effort will be 
required to come up with a good "starting point" for the 
prior model. In particular, the test frequencies and cross- 
orthogonality coefficients between analytical and test modes 
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can be used to obtain an improved prior modal model, based on 
a first-order "correction” of the original prior modal model. 
This is the "new” operation to be implemented. It will be 
discussed in more detail in Section 4. 

The procedure planned for processing response data in the 
first stage of estimation will entail sequential use of data. 
Briefly, we plan to work within selected frequency bands as 
we march along the frequency spectrum. Within each band, we 
will process sets of response data selectively. 

In the second stage of estimation where we proceed from the 
refined modal model to update the mass and stiffness matrices 
of the "dynamic model”, the critical and difficult step will 
be to define the subroatrices to be associated with the scaling 
parameters. The operation will be discussed in Section 5. 
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DERIVATION OF ESTIMATOR EQUATIONS 


This section presents the derivation of all of the equations 
which make up this Identification scheme except those which 
make up the minimum variance estimator. The derivation of 
those equations is described in References 1 and 2. The 
computer code which implements the methodology is described 
in Appendix A. ■ 

Two separate computer programs have been developed: One to 

perform the Phase I estimation of the analytic model and the 
other to perform the Phase XI estimation of the dynamic 
model. Presently the user must start each program with card 
data. This will be changed later to make the data interface 
automatic using disk or tape storage. The programs will not 
be automatically coupled, however, because the user will 
always need to evaluate the results of Phase 1 before pro* 
ceeding to Phase XX. It may, in fact, take several Phase X 
runs using different data and/or estimating different param* 
eters before an engineering* sound model is obtained. The 
methodology does not substitute for a sound understanding of 
the structure. Xt is only a tool. The better the user 
understands the structure, the more useful the tool will be 
— and visa versa -- the more he uses the tool, the better he 
will understand the behavior of the structure. 

2.1 The Minimum Variance Estimator 

The model adjustments provided by this methodology are made 
with the minimum variance estimator described in References 1 
and 2. This is a linear estimator in that it operates on the 
basis of a linearized relationship between a state vector (r) 
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and an observation vector (Y). If r_ and Y_ denote prior 

P P 

estimates of these vectors based upon analytical models , then 
the true values r and y are assumed to be related by the 
linear matrix equation 


{Y - Y 1 - ITJ {r - r } (2-1) 

P P 

If (Y> represents measured data with random error {e} and the 
♦ 

simplified notation, (y - Y^} *={y} and {r - r } » (R) is 

P P 

adopted, then equation 2-1 becomes 

{y> - (T1 {R> + {e> (2-2) 

as shown in equation 24 of Reference 2. The derivation of 
the statistical estimator equations beginning on page 187 of 
Reference 2 applies here also, except for the following dif- 
ferences 

e The analytical model is a modal representation of 
the structure, not a finite element representation. 

• The state vector may contain elements of the damp- 
ing matrix or it may contain elements of the mass 
and stiffness matrices. 

e The measured data Is the amplitude of the displace- 
ment response measured at discrete frequencies and 
locations and due to a specified set of forces* 

In essence, this formulation first trys to refine the irodal 
parameters of the analytical model while at the same time 
acting as a filter on the measured data to overcome experi- 
mental difficulties* By using frequency-response data as the 
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basis for comparison It is not necessary to excite "pure” 
modes. The use of frequency*response data also allows the 
estimation of elements of the damping matrix. 

2.2 Mathematical Models 

Before proceeding with the formulation of the minimum variance 
estimator for this parameter identification approach, it is 
necebjary to summarize the hierarchy of mathematical models 
used to represent the structural system and to define our 
notation. Four models are germane to the present discussion: 

e Finite Element Model - This model represents the 
most detailed model of the system. Created part 
by part, it may contain up to 100,000 degrees of 
freedom depending on the size, complexity, and 
importance of the structure being analyzed. 

• D ^niamic Model - This model is obtained by taking 
the finite element model, or the components of that 
model, through several stages of coordinate reduc- 
tion. This model may contain up to 750 degrees of 
freedom. At least some of these degrees of free- 
dom must relate directly to physical structural 
displacements so that a comparison to measured test 
data can be made. In the development that follows 
this model is referred to as the "u” system. 

• Analytical Model - This model Is obtained by com- 
puting a limited number of modes of the dynamic 
model and using them as the basis for an additional 
coordinate reduction. Very often for spacecraft 
structures all modes of the dynamic model (l.e., the 
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u system) below 50 Hz are selected for use in the 
analytical model. This may result in as many as 
200 and as few as 20 modes being used. The size 
(i.e.f number of degrees of freedom) in the ana- 
lytical model is the number of modes being used. 

These are the modes which should and can be verified 
with vibration test data. In the development that 
follows f this model is referred to as the system. 

e Modal Model - Prior to the Phase I estimation, the 
modal model does not exist because the coordinate 
reduction described above diagonalizes the dynamic 
model mass and stiffness matrices. There is then 
no need for another transformation. However, once 
em estimation has been performed the analytic model 
mass and stiffness matrices are no longer diagonal. 
Now a second coordinate transformation can be defined 
using the modes of the perturbed analytic model. 

The resulting coordinates are herein referred to as 
the modal model or the ”q” system* 

2.3 Equations of Motion 

Let the displacements of the dynamic model be represented by 
the vector {u(t)>. Then the equations of motion for this 
system can be represented by 

[mHu) + [CHu) + [K]{u} - {f(t)) (2-3) 


where 


[H] “ mass matrix of dynamic model 
[K] « stiffness matrix of dynamic model 


(Cl * damping matrix of dynamic model 
{f(t)> * vector of forcing functions 

For our purposes {f(t)} may be further broken down into a 
vector characterizing the force distribution {P} and a scalar 
time functionf g(t) 

{f(t)l « (Pi g(t) (2-4) 

* (P) ®g sin £2 t 

The undamped eigenvalue problem for this system may be stated 
as 

([Kl - (Ml) - 0 (2-5) 

Its solution will yield a set of eigenvalues [X] and eigen- 
vectors Using this eigenvector matrix, or a truncated 

portion of it, as a coordinate transformation, results in the 
analytical model . For our purposes, we will assume that the 
mode shapes [(^>1 have been normalized for unit generalized 
mass. We will now introduce the left-hand superscript ® to 
represent the matrices associated with the prior dynamic 
model before any changes have been implemented. Letting 

{u) « (®41 {x> (2-6) 

and premultiplying by transpose, equation 2-3 becomes 

{•41 I*41 {x> + [’4l'^{“Cj l‘»4.3 (x) 

+ 1*41 (x) - («4)^ {P> g(t) (2-7) 




•«<»v 
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but 


[•<►! 

[®4il 

[*<>1 

{P} 


Cl] diagonal 
generalized mass matrix 

[•cl 

generalized damping matrix 
[*X] diagonal 

generalized stiffness matrix 
{P} 


as a result of the modal normalization and the properties of 
the modal matrix. This results in the equations of the initial 
analytical model: 


[I] {»} +[*cl {x} + [•M {x> » {P> g(t) (2-8) 


This set of equations will now be defined as the prior analyt- 
ic model where [I] and [°X] are diagonal matrices and all 
of the matrices plus [**<{>] are a direct result of the unmodified 
dynamic model. 

The objective of the Phase I parameter estimation is to develop 
revised mass and stiffness matrices for the analytic model 
and a revised generalized damping matrix for the modal model 
which will, when incorporated into the response analysis, 
provide an improved correlation between the calculated 
responses and the measured responses. Before preceding to a 
further discussion of the parameter estimation, let us define 
the terms used in the response analysis. 


Once the analytic model diagonal mass and stiffness matrices 
have been perturbed they will no longer be diagonal. Now a 
second transformation can be defined to diagonalize these 
matrices. Let 

t'{>] * modes of perturbed analytic model 

[m] perturbed mass matrix - smalytic model 

[k] - perturbed stiffness matrix-analytic model 

and the transformation be 

{x> * I^J»] {q} 

The equations of motion now become 


where 


[11 

(q} H 

^ Ul 

{q} + 

[X] 

(ql - 


{p> 

q(t) 

(♦ 1 ® 

Im] 

m 

- (11 - 

1 * 1 ’ 

['♦ 1 ® 

[Ml 


m 

[in'* 

Icl 

m 

=• ICl - 

I'M® 

(•♦ 1 ® 

[C] 

[“ 4>1 

m 

t*]'' 

(kl 

m 

- EM - 

t*l® 

[”♦ 1 ® 

[K] 

[®«>1 



Only [I] and are diagonal. 

2.4 Frequency Response Measurements 

The displacement response of the dynamic system, denoted by 
{u}, can be represented as a transfer function between the 
response {u} and the forcing function (P}g(t). 


{uXt)} » [H(£l)]{p}g(t) 


[H(n)l also performs the same task if the problem is trans- 
formed from the time domain to the frequency domain. 

(Utn)} = {p} G{0) 


Furthermore, the transfer function matrix, for the u 

system can be expressed as a function of the eigenvectors 
(i.e., mode shapes) r the eigenvalues (X], and the damping 
matrix [c3 


or 


= [<M ($1*^ matrix notation 

( 2 - 9 ) 

m0 

“ik ’ *iq % ♦km notation 


where 


(H(n)J 


[- 


+ tcini + E^l 


] 


-1 


[z(n)]“^ (2-10) 


Derivation: (u) > [41] {x} 

= t*] [H(n>] {P}g(t) 

= [*] IH(Q)1 Ul®{P)g(t) 
=> [H(n)l {p)g(t) 


We have switched from the superscript "o" to no superscript to 
indicate that the model parameters to be used are those for 
the dynamic model as it stands at the time. Only at the 
beginning of the iterative procedure will [ 4 >] and [X] be 
identical to and [*Xl respectively. 
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Excitation frequencies Oj need not correspond to resonant 
frequencies but refer to any selected excitation frequency. 

The modulus of the complex frequency response will be used as 
a comparative basis for computed and measured response in order 
to confine the estimator to operate in real arithimetic. 

To illustrate the physical significance of the new obser- 
vation vector. Figures 1 and 2 are provided. Figure 1 shows 
a schematic representation of a free-free beam excited laterally 
at one end by a sinusoidal force g(t) » *g sin J2t, The beam 
is modeled with five transverse degrees of freedom. Three 
non-rigid-body modes and frequencies can be obtained from 
this model. Depending on the stiffness and mass distribution, 
some of these may be closely spaced. The frequency response 
characteristics over a range including the first two frequencies 
may resemble those depicted in Figure 2. Four discrete exci- 
tation frequencies are identified in Figure 2. In this case the 
effective observation vector {y} contains 20 elements including 
four vectors {u} each containing five elements. Although 
a relatively complete set of response measurements is indi- 
cated here, it is emphasized that any number of measurements 
may be considered. For example, it may be sufficient to 
include only four out of the 20 corresponding to only two 
response measurements at u^ and and two excitation fre- 
quencies and 

Intui :ively, one set of response locations and excitation 
fre^iUencies may be more suitable for calculations than another. 
Experienced judgment will be helpful in making a good selection. 
However, different sets of measurement data may be processed 
sequentially by the estimator so that one need not be restricted 
to choosing a single optimum set of data. 
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Figure 1, Finite Element Model Figure 2, Frequency Response of Model 
of a Free-Free Beaun Shown in Figure 1 
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2.5 Sensitivity Matrix Relating Observation and State 

Vectors for Phase I 

For the system identification methodology developed here, the 
observation vector is any set of discrete response measurements 
which has been chosen as data. Furthermore only the modulus 
(i.e., magnitude) of the measured responses is being used. 

The state vector, on the other hand, is composed of the 
components of the generalized mass matrix and the generalized 
stiffness matrix. Initially the generalized mass matrix is 
unity and the generalized stiffness is a diagonal matrix of the 
eigenvalues of the prior dynamic model. The sensitivity 
matrix [T] relates the state vector to the observation vector. 
The derivation of this matrix is shown below and summarized 
in Table 1. First the primary elements in the nomenclature 
need to be defined: 


Yg = luj^l * magnitude of i^^ response (2>12) 

■ i. component of observation vector 

a i, k element of lH(n)l (2-13) 


Pt. * 


Ui a «u^ + ^u. 1 

'^i ’ ^k ®ik *"k 


(2-14) 

(2-15) 


('"ik % - ®ik ^ ^k <Aik"*’k ^ ®ikX) 

(2-16) 


Yj^, the magnitude of u^ is given by 

* <“i "i> 


(2-17) 
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Equation 

No, 

(2-19) 




require 






^ requires ^ and ^ 

!iii + £i !fii\ 

^ V ) 


. 




^ik — Q -- ^ - where n defined 
Reference 6 


(2-21) 


( 2 - 22 ) 


(2-27) 


(2-37) 


’^h -\a • 51T- 


9X. 

^ - *ri *91 


(2-28) 




0 (neglect) i jt a (2-31) 

r X 9 

^i ^i ^^^i m. . i - a(2-32) 

^^rs 



Table I (Continued) 



Notes: 1) The state vector will have in__ and k _ or terms - 

X 9 79 7 9 

never both at the same time since the methodology has 
been set up to estimate separately. 
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where the * refers to the complex conjugate, the superscripts 
R and I refer to the real and imaginary components respective- 
ly, and the superscript ■* refers to the u system (i.e., dynamic 
model). 



a + b 
gm gm 


(2-18) 


The elements of [T] are given by 


Ah 




(2-19) 


where the index A- denotes the A element of the observation 
vector Y and the index h denotes the h element of the state 
vector r. Taking the partial derivative of Y using equa- 
tion 2-17; 




■ K"i “i] 


- 1/2 


3u. 3u? 

ut + u. 


( 2 - 20 ) 


but 


3u. 

±. u* 

"i 


R 


I.. 3^Ui 


u 


i 3r, 


u ^"i R Z 

As «u_. + ^u, ^ 


3r* 


i Tr 


‘i 3r, 


Therefore 


r 1 
357 * M 


- 1/2 




\x 


i ar. 


i 1 

- + u 


‘i 3r, 


(2-21) 


R I 

where and are defined in terms of the real and imaginary 
components of the transfer function and the force 

vector by equation 2-16, Thus 

- % »*ik - S 3Bik 

5^. &f. ~ ^rT“ 


3*ik , \ 3B.J, 
3r. 5rT 3r. 


Using the transformation discussed later (Section 2.7) 


or 


W » [‘^<^1 W 


^iq n ^in ^nq 


matrix notation 


indicial notation 


(2-23) 


The elements of the transfer function can be expressed as 


[HI - I"*! t*l [HI [«♦]* 


or 


E 


matrix notation 
(2-24) 


•ik - ♦in ♦nq V ‘♦ko ♦on notation 
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Thus 


— » derivative of real part of equation 2-13 
h 


]C 




in ''’^nq “qra '*'ko '^om 
n#q,m,o ^ 


* E 

n,q,m,o 


4»> «|/ 3^ 

^in ^nq qm ^ko ^om 

3r,_ 


and 


n,q,in,o _ ^ ^ 

h 


^®ik 

•— a derivative of imaginary part of equation 3 
h 


1“ 


2 ^ 1 >‘ 3'^«« H ip 

_ _ ^in nq ora ^ko ^om 

n,q,m,o ^ 


+ 2Li 3^H ®<>, tp 

^ ^ in ^nq qm ^ko ^om 

n,qfra,o ^ 


T* 


+ ^ H *A. 3tlt 

_ _ ^in ^na qm ^ko ^om 

n,qim,o ^ ^ -r- — 


(2-25) 


-13 


(2-26) 



Both emd [^] are real matrices. Only [B] is complex. 

To complete the evaluation of requires the determination 
of 3^qm and To obtain the first of these, first dif- 

ferentiate equation 2*10 and then consider the three different 
variables that r^ may be. 


3[H(n)l 


(-1) (i{a)i‘^ 




The parameters being estimated are the elements of the analytic 

model stiffness matrix, denoted by k ; the elements of the 

Ik s 

analytic model mass matrix, denoted by m ; and the elements 

t 9 

of the generalized damping matrix* denoted by E Only 
and m^^ directly influence the new model. Although 
has an indirect effect, its estimation is a separate operation. 
E „ is part of the modal model while m__ and k_„ are part of 
the analytic model. First consider 











1 31 
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0 


CASE 

then 

CASE 

then 

CASE 

then 


1) ^ * ■'rs ' 3i. 




rs 


3m.. 

®jr 

*9 


ax. 


rs 


^ -V' ♦ji *jr “ "-ri *81 

rs ^ 


( 2 - 28 ) 




and !!:ii - « 


rs 


"Sm 


jr "is 


rs 


3X. 


rs 


E 


3 > =■ 5 


*ji ^jr \s *ti 


= - *ri ^si 


( 2 - 29 ) 


h "Ts 


3k. » 3m.. 

^ = 0 and ^ ^ 


'rs 


rs 


^ - 0 ( 2 - 30 ) 

^rs 


Next consider 


(2-31) 


For r_ » m or k 

D rs rs 


H 


it 


* 0 (neglect) when ift I 


3 


3r, 




'i H 


1 / 2 . 


’ (2c when 1 


- 1/2 

Jr~ *** 


11 


(2-32) 


where 


3X. 

Is defined above (equations 2-28 to 30) 


For r. * E 
h ^rs 


3C 


11 




rs 


«ir «ls 


(2-33) 


3H 


Thus the three expressions for become 

3rj^ 


CASE 1) r. » k 

n rs 


3H \T^ *■ 

alE® - -2^ “qi On *si “im 


/ - 1/2 . \ 
ll + i) 


(2-34) 
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I 


CASE 2) r. * 

ti rs 


gin 

^“rs 


2 “qi*ri*si 


“i» h(^*° h ’'ii 


CASE 3) - C,3 


3H 


SB, 

rs 


£ 


(2-35) 


H^6. »-H H (2-36) 

qi ri 3*^ Am qr sm ' 


Finally the evaluation of ^^nq roust be provided 




♦ik 


(2-37) 


CASE 1) . k^3 




(2-38) 


3'1'i 

W 


rs 


E 


♦ik ftk 


Sj V^J - \} 


(2-39) 


* n defined in subsequent equations 
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CASE 2) r* • 

A rs 




(2-40) 




^ -s 


/ ( 1 - \ 
\ "^j ■ 


CASE 3) > 5„ 


il a 0 


ac 


(2-42) 


rs 


In equations 2-16# 21, 25, 26, 33 to 35, 39 to 42, we now 
have all that is needed to calculate the elements of the 
sensitivity matrix, 
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2.6 


Modification of the Analytic Model (Phase I) 


The prior analytic model, as represented by equation 2-8 has 
the mathematical characteristics sho%m in equation 2-7, This 
model, however, does not provide an exact match of the 
measured test data. If it did, there would be no need to 
improve it. Therefore there must exist other sets of modes 
[0] and frequencies iX] which will provide a better repre- 
sentation. Therefore our goal is to develop a new analytic 
model mass matrix, [m] ; a new analytic model stiffness matrix, 
tk} ; and a new generalized damping matrix, [CK 

These new matrices will not, however, satisfy the require- 
ments of equation 2-7: 

[♦1® I'M) [*1 I« til 
[♦l'' t'Kl t*I * [XI 

The approach to be used to estimate these new matrices is as 
follows: 

A) Calculate the response of the system at the measurement 
locations using equation 2-8. [H(n)l is given by equation 
2-10 and [01, [X], and [^] are taken from the unperturbed 
dynamic model. 

B) Calculate the sensitivity matrix [T] and the "observation" 
vector. The effective "observation" vector is the dif- 
ference between the measured responses and the calculated 
responses (equation 2-2). 

C) Use the minimum variance estimator to provide new mass 
and stiffness matrices for the analytic system. Input 
to the estimator consists of the "observation" vector, 


- 24 - 


) 


r 






the sensitivity matrix, and the previous estimate of the 
mass and stiffness matrices. The new equations of motion 
for the "x" system are 

[mHx} + Ie]{x) ^ (JcHx} - {P)g(tJ 

D) Solve for the eigenvalues jX] and eigenvectors [t(»] of the 
modified "x” system. Normalize these modes such that they 
have the characteristics: 

Im) * 1 13 diagonal 

[^1^ Ik] [(}»] » IX] diagonal 

E) The revised eigenvectors for the dynamic model, the u 
system, are given by 

[♦1 - I®*] m 

F) Revise the damping matrix [?] (from equation 2-10) to 
reflect the new eigenvalues while retaining the same 
damping ratios as were assumed for the initial dynamic 
modelf 

G) Calculate the response using equation 2-8 but with the new 
eigenvalues [X] and new dynamic model modes Ul. 

H) Repeat the above steps until convergence is obtained. 

I) Lastly, perform a similar iterative scheme to estimate 
elements of the damping matrix, 


*It is assumed here that damping parameters are 
estimated by a separate operation (see Step 1) . 
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2.7 Sensitivity Matrix Relating Observations and State 

Vectors for Phase II 

The next step in the two-phase system identification procedure being 
developed is to generate improved mass [H] and stiffness [K] 
matrices for the dynamic model, the u system. The procedure 
involved here is similar to that of Phase I, if one con- 
siders the "prior” model now to be the dynamic model instead 
of the analytic model and the "test” data to be the refined 
parameters of the analytic model instead of the response 
data. 

It is impossible at this point, however, to estimate all of 
the elements of the [M] and [K] matrices since these matrices 
may be as large as 750 by 750. Instead these matrices will 
be broken down into a number of sub-matrices each of which 
can be scaled individually. 

The formulation of Phase II of the system identification 
procedure is shown below. It envisions that the variable 
parameters of the dynamic model mass and stiffness matrices 
will be scaling constants such that 

[®M] =* dynamic model mass matrix 
n 

* [“Ml ®h (2-43) 

h=l 
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and 


[•K] • dynamic model atiffnesa matrix 

^ l«Kl 
h«n+l 


C2’-441 


where 


I*H3 ^ Portion of mass matrix which, does not vary, 

[*K] « Portion of stiffness matrix which does not vary, 
* Localized portion of M which is to be scaled, 
[**Klh * Localized portion of K which is to be scaled, 

= Scaling parameter for either [M]j^ or [Klj^ and 
N a Total n\mber of variable degrees-of-freedom in 
[M] and [K] combined 


In indicial notation equations 2-43 and 2-44 become 

n 


°“rs=‘®rs*E “h'“rsh 
h=l 


'«rs ^ 


.N 


E “h"rsh 
h=n+l 


(2-43) 


(2-44) 


In the Phase II "prior model", it will be assumed that = 1 
for all h.- Variances will be specified independly for each 
Along with the revised parameter estimates from the 
analytical model will be a corresponding covariance matrix, 
containing the variances and covariances of the Phase II 
"test data". Thus it follows that for the Phase II parameter 
estimation, the state vector will be 
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{r} - {«} 


while the first portion of 
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At the completion of the Phase I estimation, a non-diagonal 
generalized mass matrix [m] and a non-diagonal generalized 
stiffness matrix [k] have been obtained. The objective of 
the Phase II estimator is to use these two matrices as test 
data in developing improved dynamic mass and stiffness 
matrices, [M] and [K]. "In addition, an improved set of mode 
shapes [i^l have resulted from Phase I. These modes can be 
used to develop a relationship from which the Phase II 
sensitivity matrix [T] may be obtained. 


We have defined 


[(f) } a improved modal matrix of dynamic model 
® prior modal matrix of dynamic model 
[ip] 3 modal matrix of revised auialytic system 
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Thus 


M imi m * m 
[k] m * IM 


(2-47a) 

(2-47b) 


What we want to obtain are revised [M] and [R] matrices such 
that 

Ul* IM] Ul = [X] (2-48a) 

[* 3 ® [K[ I4J - [I] (2-48b) 

We also have the relationship between [(^] and [ifr] 

W - ["(frl m (2-49) 

First substitute equation 2-49 into equation 2-48b 

tn)* t’41® tkl ['♦! I'M - IX) (2-50) 

From a comparison o£ this relationship with equation 2-47br 
we obtain the relationship 

tk) - CKl rol (2-51) 

where [Kl is the presently unknown matrix which we wish to 
generate. Substituting equation 2-44 into the above relation- 
ship provides one o£ the needed relationships. 


[kl - 


y h=n+l f 




(2-52) 


or in indicial notation 


^ij 


■D V*rshr*sj 

r s \ h«n+l / 


Using equations 2-47a, 48a, and 49 and a similar argiiment 
the relationship for [m] may be obtained: 

n 

n r OMl 

h 


[m] - 


l*M] 

h=*l 


[*4»1 


(2-53) 
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or 


m 


ij 


E •♦rlKs-E V«rshV*sj » 

r ,3 \ h*l / 


-53) 


The elements of the sensitivity matrix IT] for Phase ;i param- 
eter estimation are still of the form given by equation 2-19: 




CASE 1 


BY I 
30t^ 


“ij 

E 

r,8 


r,3 


:‘.i( £”•«„.?%) 


(t (|) 


sj 


X) ”rsh “♦sj for i - 1 to n (2- 


54) 


BYi 


0 for 4 » n+1 to N 


CASE 2 


“"ij 


S-E 


r,s 


\ t=n+l 


E '♦rl‘*rsh •♦aj ‘ ** 


(2-55) 


r,s 


3 a. 


0 for i » 1 to n 
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Nocice that the sensitivity matrix is not a function of the 
scaling parameters, ot^. Only the prior modal matrix 
and the unperturbed dynamic model mass [M] and stiffness IK] 
matrices enter into its calculation. This, of course, is 
a direct consequence of the linear relationships between K 
and k, and between°M and m. Because of this linear relation- 
ship, the iterative procedure will converge in only one step. 

2.8 Modification of the Dynamic Model (Phase II ) 

The goal of the Phase II estimator is to develop a set of 
scaling parameters, a, which will improve the mass and stiff- 
ness matrices of the dynamic model. It still remains the 
responsibility of the analyst to select a set of submatrices 
which when multiplied by scaling factors will improve the model. The 
methodology described heretofore does not, unfortunately, 
provide any insight as to how these .submatrices should be 
chosen. The analyst must select subma trices which he thinks 
might be successful or enlightening based on his experience 
and knowledge of the structure being investigated. In general, 
a number of trial configurations may be run before a useful 
and realistic modified model is obtained. 

The basic approach used to estimate the new dynamic model 
mass and stiffness matrices is as follows; 

A) Complete Phase I to provide ** observation" data in the 
form of generalized mass [m] and stiffness [k] matrices 
for the dynamic model and a covariance matrix for these 
elements. It must be emphasized that the generalized 
mass ^uld stiffness matrices for the dynamic model (the 
u system) are the mass and stiffness matrices of the 
analytic model (the x system) 
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o) Calculate the sensitivity matrix [T] . This need only be 
done once because the sensitivity matrix is only a 
function of the unperturbed prior model. 

C) Determine the "calculated*' generalized mass and stiffness 

matrices using the prior modal matrix and the unper- 

turbed dynamic model mass and stiffness matrices. 

D) Form the "observation" vector as the difference between 
the elements of the "calculated" [m] and [k] matrices and 
the "observation" {m] and [k] matrices. 

E) Use the minimum variance estimator with the observation 
vector# the sensitivity matrix and the prior estimate 
of the scaling parameters to provide a new estimate of 
the scaling parameters. The prior estimate for each 
scaling parameter is 1.0. 
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EXAMPLE PROBLEMS 


At the start of the identification procedure 

[ml - (•♦1* (“Ml [•♦! = [*111] = [II (3-1) 

[!0 • [•♦l’^ ['Kl [»♦] - [*Icl = [Mv.1 (3-2) 

where [^4>] * modes of prior model 

{®M] - mass matrix of prior model 
[*K] = stiffness matrix of prior model 

Ideally, for the analytical examples present in this section, 

[m] and [kj should converge to 

[ml - ["iH® [Mgl ["itl (3-3) 

[kl - [‘♦l'^ [%1 [‘*1 (3-4) 

where [M 1 and [K ] are the known (in the real world they 
e e 

will be unknown) mass and stiffness matrices of the correct 
model. Phase X of the procedure attempts to find the correct 
[m] and [k] matrices while Phase II uses the results of Phase I 
to identify [Ml euid [K 1 . Only the upper triangular elements 

xS c 

of [m] and [k] are sought (i.e., symmetry is assumed). 

This section discusses the four example problems investigated 
to date. Before the individual problems are discussed, a 
couple of general comments should be noted: 

• In the following discussion, "converged® means that 
the responses calculated with the modified [m] and 
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[Ic] arts within the specified tole. ance band of the 
exact values. 

e **Fit improved" means that [m] and Ik} converged to 
values which significantly improved the prior model 
but did not fall within the convergence tolerance 
of desired values. 

e Except for Run 3-3, the "test data" used were exact 
values obtained from an analytic solution of the 
exact model. For Run 3*3, random errors with a 
standard deviation of S% were introduced into the 
"Test data". 

• A variance, or standard deviation, is specified for 
each element of [m] and tic] being estimated and for 
each observation. 

• On four of the runs which did not converge, only a 
portion of the [m] and [k] matrices were being 
estimated. Neglecting the farthestmost off-diagonal 
term created difficulty for the 3 degree-of- freedom 
problem. 

e The terms "exact", "test", and "observation" all 
refer to the same set of parameters. The "exact" 
model is what we wish to eventually ^u^rive at. It 
is the model used to artifically generate "test" 
data from which "observation" points were selected. 

The basic conclusion which we draw from these runs is that the 
estimation methodology works well provided the "prior" model 
is sufficiently close to "exact" model. 
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3.1 


One Deqree^of*Freedom Problem 

This, the simplist of all possible problems, has the great 
advantage that simple relationships can be derived for the 
sensitivity matrix and a hand solution can be obtained. There 
is, however, no Phase II estimation required for this problem 
because = 1* The equation of motion is 

Mu + Cu ■»* Ku » sinSlt u system (3*5) 

u + 2CW u + Xu » 5 ^ sinnt x system = u system (3-6) 

O M 


The transfer function is: 

H(!J) » — = A + Bi (3-7) 

(u ^-0^) + 

O 0 

2 2*^ 

A « (w 

o 

where 

X =« Wq 

The sensitivity matrix is given by: 


3Y, 

» 3 i 

Ih 3r. 


[uu*] 


- 1/2 


^U. a^U 

^fh '■^h. 


(3-8) 


Fa 3A ^ B 3B 1 





where H » tA^ + - /m 

p = ^ + * ^ .+ 0£ * ^ 

First Parameter: ■ k 

_ 

3k ■* 



(3-10) 




The operations o£ the computer program to calculate the sensi- 
tivity matrix elements were verified with this example problem. 
When the exact model and the prior model were assigned the 
following values: 


• Exact Model; K=*1.0 X^l.O 

M » 1.0 4> =■ 1.0 

5 » 0.1 
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• Prior Model 



the estimator converged to within 0.1% of the desired values 
in 6 iterations. 



3.2 


Two Degree-of-Freedom Problem 



The two degree-of- freedom example 
is represented by the adjacent 
sketch. The artifically generated 
"test” data is based on the 
following parameters; 

kj^ - 1.00 kj = .0125 

* 1.00 m2 * .20 

» .10 ^2 ■ *05 


The resulting "exact" frequencies and modes are 


» 0.1952 rad/ sec 
*^>2 * 1.0042 

W *r .02938 .99957' 

L2. 23510 -.06570 


The prior model used to investigate the behavior of the Phase 1 
estimator had the following characteristics; 


kj^ =« 1.00 

ra^ * 1.0 
- .10 
<a^ • .1992 


k 2 * .008 

IB 2 * 0.2 
C 2 * .05 

0)2 * 1.0042 


In all cases we restricted the investigation to constant mass 
and constant damping* The only difference between the "prior” 
model and the "exact” model was in the stiffness parameters. 
Eight runs were made. These are summarized belows 





Risulta 


Run llo. TfSt Ringe of Mo. Tt»t ® 

No. Frtq. Test Fr«q. Pofnts ptr l!cj M 

(rid/tic) Frtquency 


t.2 

6 

.19-1.06 

2 

.C02 

.000001 

Converged to 




.002 

.000001 

within 5X In 





.100 

.000001 

4 Iterations. 

3 

6 

.19-1.06 

at 

<i 

.002 

.002 

Old not Improve 





.002 

.002 

fit In S 

* 




.100 

. 1/02 

Iterations. 

4 

10 

.10-1.0 

1 

.002 

.002 

Converged to 





.002 

.002 

within Zi In 





.100 

.002 

4 Iterations. 

5 

10 

.10-1.0 

1 

.002 

.002 

Only estimated diagonal 





.ICO 

.002 

elements of [m] and [k] 
Estimator diverged. 

6 

10 

.10-1.0 

1 

.002 

.00001 

Converged to within 





.002 

.00001 

2X In 4 Iterations. 





.100 

.00001 


7 

10 

.10-1.0 

1 

.002 

.00001 

Eld not converge to 





.ICO 

.00001 

desired model but did 
Improve fit of data. 

8 

c 

.19-1 .C6 

2 

.002 

.00001 

Only estliiiated diagonal 





.100 

.COCCI 

elements of and [k] 

Fit Improved, 

* Here 0 ^ 

= variance 

. Sequence Is 1- 

1, 1-2. 

C-Z. 

- 



Several typical frequency-response plots are shown in Figures 
3, 4, 5, and 6. 

On Run 1,2 the estimator was provided with six test points at 
each node. The tolerance on [k] was set very wide and the 
tolerance on [m] was set very tight. In 4 iterations the rms 
error between the calculated and the test responses was re- 
duced from 1.29 to 0.05 with all of the calculated responses 
within 5 percent of the "test" data. (Figures 3 and 4). 

Run 3 is the same as Run 2 except that the prior estimate of 
the tolerance on [m] was set very wide. As a result the esti- 
mator overshot the desired values and then had trouble finding 
its way back. After 10 iterations it was not able to improve 
the fit between the calculated and the test response. 
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Figure 4. Frequency Response Plots, Two Degree -of -Freedom 
Example, Mode 2 
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On Hun 4 the estimator was only given observation data on one 
node, node 2. However, the tolerance on both [m] and [k] was 
left very loose. The estimator was able to match the ^test” 
data to within 2% after only 4 iterations* The Im] and [k] 
were not as good as on R\m 2 however. The new natural freguen- 
cies were almost identical to the desired values but the new 
mode shapes were not as good. Because data was provided for 
only one node, the new model matched that data very well, but 
the response of the other node had drifted off. This demon** 
strates well the fact that a set of parameters which fits the 
test data well but are unrealistic may be obtained if the test 
data is not representative of the entire model. {Figures 4 and 5). 

On Run 5 we tried to estimate only the diagonal elements of 
[m] and [kl with wide tolerances on both the prior estimate 
of {mj and prior estimate of [k] . In addition, data was pro** 
vided only for node 2. These conditions were just too severe. 

With no coupling between modes and no observation data on node 1, 
the estimator has no way to estimate the needed mode shapes. 

After 5 iterations, one of the diagonal elements of the mass 
matrix went negative. 

Run 6 is identical to Run 2 except that observation data was 
provided at only one node, node 2. The estimator was able to 
generate a set of parameters that matched the given response 
data to within 2% after only 4 iterations. But just as on 
Run 4, the new parameters were not as accurate because the 
response on the other node was not controlled. 

Run 7 is similar to Run 6 with a similar conclusion. After 10 
iterations the rms error between observation and calculation 
was reduced from 1.47 to .02. But because one node was 


uncontrolled and because the off«diagonal terms were not 
estimated# the new model was not realistic even though the 
given test data was well matched (Figures 5 and 6). 

Run 8 is similar to Rton 7 with only the diagonal elements being 
estimated but with data provided at both nodes. The estimator 
was not able to match the test data as well as on Run 7 or the 
desired model as well as on Run 2. But a much more realistic 
model was obtained. The fact that data was provided at both 
nodes prevented the estimator from pulling away from the desired 
parameter but the lack of coupling (only diagonal elements 
were estimated) prevented obtaining a highly accurate set of 
parameters (Figures 5 and 6). 

For execution of Phase II of the parameter estimation# the 2 
by 2 matrices were so partitioned that each element was 
separately estimated: 


Prior Models 





1.0 0 
. 0 0. 

" 4 : 

IKJ .f® 

“l+Oj 


lo 

oj 

. 0 OJ 1 

Exact Model: 



>«el 

.0 0 
0 0.2 

j 



.0125 

.0125 

-.0125 

.0125. 



008 -.008 
008 .008 
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Figure 6. Frequency. Response Plots, Two Degree of Freedom 
Example, Node 2 
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• Desired Generalized Stiffness Matrix [K] 


Ikl - 


.0617917 

-.0103875 


-.0103875 
1.01321 . 


• Prior Modes; 


.018470442 

.2.23568652 


.999829407] 
-.04130120 J 


• Desired Generalized Mass Matrix 




[M] l.iH: 


[m] 


i.o 0 ] 
. 0 i.oJ 


• Desired Scaling Parameters 

» 1.0 

02 * 1.0 

s 1.004464 
0^ * 1.562500 


Phase I results from Runs 6, 1 , and 8 were all used with the 
Phase 11 estimator, in all cases, the estimator produced a 
set of scaling parameters such that 

[•*]'* [Ml I"*! - [B] 

[•♦I* IKJ [•♦! - [kl 

converged within the specified tolerance. The closeness of the 
new [M] and [K] matrices to the desired tM.] and [K 1 matrices 

G G 

depended, of course, on the fidelity with which [m] and [k] had 
been estimated. Of course, if fewer scaling parameters had been 
used, an exact match would not have been possible. 
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3.3 


Three Degree*of Freedom-Problem 



The three degree-of-freedom 
example la represented by the 
adjacent sketch. The artificially 
generated "test* data is based on 
the following parameters: 


= 47.0, 


31.0, 

*^3 " 

39.4784 

» 0 . 8 , 

”2 “ 

1 . 2 , 

m 3 . 

1.0 

» 0.5, 

«2 * 

.05, 


.05 


The responses obtained with these parameters are called **test 
data "observation" data, or the "exact" model to indicate the 
various ways this methodology can be approached. The exact 
natural frequencies are: 


- 2.700 rad/sec 
0^2 = 8.284 
0 )^ = 10.945 

and the exact modes: 


[01 


.25436 

.59217 

.72625 


.63390 

.47362 

.63981 


-.88481 

.50913 

-.25025 


The prior model upon which all of the example -three runs began 
had the following parameters: 
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t 


I 


r 


'‘l- 

39.5 

k 2 « 39.5 

ks « 

39,5 

■ 

1.0 

m 2 * 1*0 

m 3 - 

1.0 


.05 

^^2 - -05 

^3 " 

.05 


This yields system natural frequencies of 2,?96r 7.835, and 
11.322 rad/sec respectively. In this example we have both the 
mass and the stiffness of the prior model incorrect. Seven 
runs were made. These are summarized below. 


Run 

No. 

NO. Test 
Freq. 

Range of 
Test Freq. 
(red/ sec) 

No. Test 
Points per 
Frequency 

o^for' 

CO 

c»fbr* 

M 

Results 

1 

9 

2.5-11.3 

3 

6* 10.0 

.10, 5*.01 

Converged to within 
51 In 9 Iterations 

2 

9 

2.5-11.3 

3 

5* 10. 

.10, 4*.01 

1-3 tern not estimated, 
estimator diverged. 

1 

9 

2,5-11.3 

3 

5* 10.0 

.10, 5*.01 

Fit Improved. 

4 

3 

see discussion 

3 

5* 10.0 

.10, 4*.01 

Estimator diverged 

5 

3 

see discussion 

3 

6* 10.0 

.10. 5*.01 

Estimator diverged. 

So 

3 

see discussion 

3 

6* 10.0 

.10, 5*.01 

Fit improved. 

Sb 

3 

sue discussion 

3 

6* 10.0 

.10, S*.01 

Fit Improved. 


# him 0 ^ = vtrluict. $«<)uenct ts 1-1, 1-2, 1-3, 2-2« 2-3, »nd 3-3. N* aeans rcptat N tines. 


On Run 1 the estimator was provided with 9 test points at all 
three modes. These test points bracketed each exact/prior 
natural frequency. The tolerance on both (ml and [k] was set 
wide. In 9 iterations, the rms error between calculated and 
exact response was reduced from .040 to .002 and a good approxi- 
mation of the desired [m] and [k) matrices was obtained. 
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Run 2 is identical to Run 1 except only the diagonal and first 
off-diagonal terms of [m] and [k] were estimated* Just as 
happened in a similar situation on example two, the estimator 
overshot on its first iteration so far that it was not able 
to find its way back. 

Run 3 is identical to Run 1 except for a new wrinkle; random 
error was introduced into the test data. In this instance, 
the estimator was not able to bring all of the calculated 
responses within 5% of the test values but it was able to con- 
verge on to a set of parameters which significantly improved 
the comparison. The rms error was reduced from .040 to .007. 

Run 4 presented a very severe problem to the estimator because 
only three observation frequencies were used: one at each of 

the natural frequencies of the exact (i.e., the test) model. 

The sensitivity matrix IT] is closely related to the transfer 
function matrix [H] . When the observations do not bracket the 
natural frequencies of the prior model, the sensitivity of the 
estimator is somewhat unpredictable. In this case, the estima- 
tor overshot so far that the diagonal elements of the estimated 
matrices went negative. The number and location of the obser- 
vation points are very important in determining whether or not 
a successful run will be accomplished. Run 5 was made to see 
whether addition of the 1-3 term to the estimation would have 
any beneficial effect. It didn’t. 

Following Runs 4 and 5 of this problem, a two-way step retarder 
was incorporated into the methodology. This retarder checked 
whether or not the new estimate had a negative on the diagonal, 
or an increased rms error. If either was found, it halved the 
step and tried again. Run 5 was then repeated as Run 5a. The 
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lack of data still presented a serious obstacle. But after 
20 iterations or partial steps, the rms error was reduced 
from .21 to .011. Then the retarder was changed from a 1/2 
to a 1/4 step and Riin 5 repeated again. Again progress was 
slow but positive. After 27 steps the rms error was reduced 
from .021 to .006. 


The [m] and [kl matrices generated on Run 1 were used as in- 
put to the Phase II estimator. Sufficient scaling parameters 
were used to provide the opportunity for an exact match. A 
close fit was obtained after only a couple of iterations 
(see below) • 


Prior Models 

[m] 


0 
0 
0 

1.0 


[kl 


0 0 1 

0 0 1 + 
0 l.oj 



+a. 


0 0 0 

0 1.0 0 

0 0 0 




+a. 


+a. 


'39.4784176 -39.4784176 0 

-39.4784176 39.4784176 0 

0 0 0 


+a, 


0 

0 

0 


39.4784176 -39.4784176 

-39.4784176 39.4784176 


Exact Model 




0.8 0 0 

0 1.2 0 

0 0 1.0 


m 

[K^] = 78.0 -31.0 

70.47841761 

SYMMETRIC 

0 

-39.47841761 

39.47841761 

Scaling 

Parameter 

Desired 

Obtained 


0.8000 

0.8017 


1.2000 

1.1809 


1.1905 

1.2039 

^4 

.7852 

.7664 


1.0000 

.9934 
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3.4 


Two Deqree^of-Freedom Problem; Close Modes 



To investigate the effects of 
close modes on the estimator 
performance, the two degree-of« 
freedom model was adjusted to yield 
two resonant frequencies within 20% 
of each other (Figures 7 and 8) . 

The artificially generated "test" 
data is based on the following 


parameters: 



» 33.333 

^2 - 

1.0 

« 28.83 

* 

.8649 

5“^ * . 05 

^2 • 

.10 


The resulting frequencies and modes ares 


=» .98617 rad/sec 
Wj “ 1.17214 


t<|i] =*r.l2588 

.13726 “1 


[.79246 -.726767J 


The sinusoidal forcing 
magnitude of 1 pound. 

function was 

applied at node 1 with a 

The prior model used here had these 

characteristics : 

* 86.715 

^2 * 1.00 


3 75.00 

m2 * .8649 


« C, - .05 

n 1 

«2 - -lO 
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Figure 7. Frequency Response Plots, Example 4 - 
Close Modes, Node 1 


- 54 - 





FREQUENCY (RAD/ SEC) 


Figure 8 
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Frequency Response Plots, Example 4 - 
Close Modes, Node 2 




E 


which yielded natural frequencies of 1.0191 and 1.1346 rad/sec. 
Three runs were made: ^ 


Run 

No. 

No. Test 
Frtq. 

Range of 
Test Freq. 
(rid/sec) 

No. Test 
Points per 
Frequency 

0^ for 

2 

o fbr 

M 

Results 

1 

10 

.9-1.25 

2 

all .0025 

all .0025 

Unsuccessful - no 
step retarder. 

2 

10 

.9-1.25 

2 

all .0025 

•11 .0025 

Unsuccessful - one-way 
step retarder. 

3 

10 

.9-1.25 

2 

all .0025 

•11 .0025 

Set discussion. 


The first run was performed without the step retarder discussed 
in the previous example. As a result the second iteration 
produced such a large change in Cm] that a negative diagonal 
element was generated. Next a one-way step retarder was added 
to prevent large increases in the rms error between the calcu- 
lated and the measured (i.e,# test or exact) responses. This 
enabled the estimator to have a better first iteration but 
one which still overshot the mark. Thus it still produced 
a negative diagonal mass element when it tryed to return 
(Run 2} . Next the two-way step retarder was added. This 
prevented the estimator from "blowing-up" but it did not en- 
able it to effectively improve the model. Apparently, the 
prior model is so far off in its mass and stiffness that the 
sensitivity matrix elements are not representative of the true 
model. This is true even though the frequencies are well 
matched. 

Phase II of the estimation procedure was not performed on 
this example. 
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4. 


METHOD FOR IMPROVING PRIOR ANALYTIC MODEL 


One of the difficulties pointed up in the foregoing examples 
is that when the prior analytic model is not sufficiently 
close to the physical system, in the sense that computed 
frec[uency response does not match experimentally measured 
frequency response well enough, then the estimation procedure 
may not converge properly. Another difficulty pointed up 
is that of identifying which of the off-diagonal elements 
of [m] and [k] need to be included in the parameter vector 
being estimated, i.e. which elements of these two matrices 
need to be adjusted to achieve a proper fit of the model 
to the data. This section documents an analysis aimed at 
resolving these difficulties. 

4.1 First Order Correction 


Two observations may be made at the outset: (1) there are 

certain measurement data available from modal survey tests 
which are not explicitly used in the estimation procedure 
described in Section 2, namely mode shapes and frequencies; 
and (2) previous experience with MOUSE [1] indicates that 
when such information is used in a linear estimator to revise 
a prior model, most of the change to the parameters of the 
prior model is made in the first step of the iterative 
estimation procedure. These observationi^ suggest that 
direct use of measured modal data may provide the key to 
conditioning (or improving) the prior analytical model so 
as to enable it to converge to the measured response data 
using the methods of Section 2. The implementation of some 
procedure to affect this conditioning, in conjunction with 
the existing ccxnputer programs, has the potential capability 
of combining the distinct advantages of two different approaches. 
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The first approach utilizes explicit modal data (where 
available) as a basis for improving the analytical model; 
the second approach utilizes sinusoidal response data where 
"pure modes" are not required. The primary advantage of the 
first approach is that where good modal resolution has been 
achieved, numercial computation appears to be more stable. 
However, only mass and stiffness parameters can be esti- 
mated, no damping. The advantages of the second approach 
are (1) damping parameters can be estimated as well as mass 
and stiffness, and (2) "pure modes" are not required. Thus, 
in the case of closely spaced modes where experimental iso- 
lation of the modes may be hard to achieve, sinusoidal 
response data reflecting multi-mode participation can be used. 

Let us assume that our "prior dynamic model" is given in 
terms of the mass and stiffness matrices [*’M] and'[^K]. 

Here again, the left-hand superscript notation is used to 
avoid double subscripting later on. The corresponding 
eigenvalues, and eigenvectors, are obtained by 

solving the eigenproblem 

([®K] - » {0} (4-1) 

Now suppose that there exists a "true dynamic model" in the 
sense that it produces a set of eigenvalues, Xj, and eigen- 
vectors, which agree with measured test data. Let 

that model be represented by the mass and stiffness matrices 
[M] and [K] respectively. Thus 

([K] - XjtMl){<|>}j - {0} (4-2) 

We can relate the parameters and modal characteristics of 
the prior model to those of the "true" model by 
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[K1 

- l«Kl + [AK] 

(4-3a) 

[Ml 

- [«M] + [AMI 

(4-3b) 

Xj 

» *X^ + AXj 

(4-30 


- {A(fr>^ 

(4-3d) 


Then, substitution of equations 4-3 into equation 4-2 gives 

[(“K + AK) - + AX^) («M + AM)j{®4>j + A(J>j} » {0} (4-4) 

Subtraction of equation 4-1 from 4-4 leads to the first-order 
approx ima t ion : 

[®K - »Xj«M]{A(|>}j + [AK - ®XjAM - AX^MH®<tt}^ * (O) (4-5) 

Following the same procedures used in deriving eigenvalue 
and eigenvector derivatives 13], we first premultiply 
equation 4-5 by When i = j, we find that 

AXj - {*c^}^tAK - «XjAMH®4>}j (4-6) 

when i j, we find that 

^[*K - “Xj®MHA<(>}j + lAK - ®X^AM) {0} 

The key step at this point is to recognize that {A^} may be 
expressed as a linear combination of the prior model mode 
shapes Thus: 

{A*fr}^ 2[«4>KAn>^ (4-7) 

Substitution of this form into the previous equation results 
in the following for i ^ j : 
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(4-8) 


Note the similarity in form between equation 4-6 and equation 
4-8. 


One additional piece of information is required to complete 
the picture, so to speak, i.e. we need a corresponding expres- 
sion for An^j, the diagonal elements of the matrix [AtiJ. Clearly 
equation 4-8 is meaningless whenever i =» j. The additional 
piece of information is provided by the normalization condition 
imposed on the eigenvectors Let us require that 

{d)}j s 1 (4-9) 

Since we have assumed originally that 

we find that substitution of equation 4-3 and 4-7 into 4-9 
yields tha final equation 

Arijj - [AM] {•((►> j (4-10) 

At this time we recall that our objective is to improve cur 
"prior analytic model, "originally specified in terms of reduced 
mass and stiffness matrices [°m] and [°k] * Let us define 
a first-order correction to these matrices in terms of the 
incremental adjustments to its elements 

s Am^^ (4-lla) 

{®<!>>JtAK]{"<fr}j i Ak^j (4-llb) 
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We further recognize the symmetric form of [m] and [k] in 
general, so that 






Thus, if we are to determine an improved prior analytic 
model in terms of say [m] and [k] , we find that we may use 

[m] - [®m] + [Ami - [I] + [Am) (4-12a) 

[kl - [’k] + [Akl - + [AkJ (4-12b) 


In general, if [ra] and [k] are of order n x n, we will require 

(n^ + n)/2 terms to define each of [Am] and [Ak] , or a total 
2 

of (n + n) terms for both of them. With reference to equa- 
tions 4-6 and 4-8, and 4-10, we find that if we know An^^^ for 
al3 i and j, each ranging from one to n# and in addition knew 
AX^ for all j, then we can solve these three equations directly 

for Am, . and Ak , . . 

i3 


To complete the derivation, we recall from equation 2-23 that 
and ^recognize from equation 4-7 that 


(Anl 5 w - m 


(4-13) 


Finally, we recognize that the cross-orthogonality between 
the measured "test modes" and the original modes 

may be expressed in terms of 

[Jj - - iJl (4-14) 
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Therefore, 


= '"ij - hi 


(4-15) 


where 6., is the Kronecher delta, 
ij 

A 

In summary, given the test frequencies, we may compute 
"test eigenvalues" Xa * u. , Given the test mode shapes, 

A J 3 A 

we may compute the cross-orthogonality coefficients i|>. . 

j ^ j 

From this information we can then compute 


AX . 

1 


X. - "X. 
j 3 


Using equations 4-6, 4-8 and 4-10 we may precede to compute 
Am. . and Ak. .. Then from 4-12 we obtain the mass and stiff- 

^3 Ij /I, A 

ness matrices, tm) and [k] , for an improved prior analytical 
model. 

The computations involving equations 4-6, 4-8 and 4-10 
require only the solution of pairs of linear algebraic 
equations. To obtain the diagonal terms, equations 4-6 
and 4-10 are used, recognizing that 




j 






Thus 4-6 and 4-10 become 


AXj =* Ak^j - ®XjAm^^ 


^^ij * ’i'ij " ^ * I^"*ij 


(4-16a) 

{4-16b) 
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In a similar manner, equation 4*8 yields 
“ {^^tj " 

4.2 Criteria for Banding [ml and [k 

It may be recalled that the original (priory "dynamic model” 
was defined by t®M] and L®K1 with corresponding equations of 
motion 

[®Ml{u> + [®C]{u} + l»Kl{u} * {f(ty} 

Transformation to a modal coordinate system was defined 
such that {u> * resulting in what we have called 

the prior "analytic model." Equations of motion in the 
x-coordinate system are 

[®m]{x} + [®cHx} + (®k){x} - t®(|»]'*^{f(t)} 

We have denoted the (unknown) "true" values of mass, damping 
and stiffness matrices in these two coordinate systsms as 

[M] , [C] , [K] Dynamic Model 
[ml , [c] , [kl Analytic Model 

where 

S lol 

f+i’icir*! s tc] 


(4- 17a) 
(4-17b) 
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We have also derived in Section 4.1^ approximate values 

* At 

for [m] and [k] which are denoted by [m] and [k] respectively. 

Let us now consider the following equations of motion 

tioUx) + [SHx) + [kj{x) - 

- [•4ij'{P}g{t) (4-18) 

We have not yet defined Cc] # but introduce the notation for 
convenience here. Transformation to the frequency domain gives 

([k] - Ji^[m] + * [«<!>]’^{p}G(in) = {p}G(i«) 

' ' (4-19) 

It is furthermore convenient to introduce the notation 

X(iSi) = X(iO)/G(ifl) (4-20) 


NOW let us recall that 

Era] =*» [l] + [Am] 

[kl » + [Ak] 

and for the sake of consistancy, let 
[Z] » f’c] + [Ac] 
where [**0 3 is diagonal such that 





(4-21a) 

(4-21b) 

(4-21C) 


and is the prior estiamte of the critical damping ratio 
for the mode. 


Zt follows that we may define the dynamic impedance matrix 
in the x*coordinate system as 
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IZ<iQ)J * [k] - a^(m] + mtcl 


<4-22) 


where the tilde over the Z is is used to denote the fact 
that elements of ? correspond to the x-coordinate system,* 

Recognizing the form of equations 4-21, we may express 
as the sum of diagonal and nondiagonal complex 
matrices. 

[Z(itl)1 » ['Zdn)] + [AZ(in)l (4-23) 

where 

[®Z(i«)] * + i£lf*c] 

[AZ(ift))- [Akl - + ifltAc] 

We can then use |[®Z(ti3)3 to generate a complex scaling trans- 
formation , 

[d(1o)] s (4-24) 

such that a new coordinate system, y, is defined by (x) =■ 
[dJ{y>. Then equation 4-19 becomes** 

{[I] + [0][45][D])jr| - (4-25) 

The purpose of making this scaling transformation is, that 
now the numerical significance of Am^^ and Ak^^ may be 
examined, by evaluating the corresponding term of the matrix 


*To be entirely consistent in notation, the carat ('^} should 
appear over the 2, but is omitted for convenience. 

**The T" notation in equation 4-25 is used in the same sense 
as the in equation 4-20* 
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[d][A 2 ][d] in equation 4-25. If any term of this matrix 
is small compared to unity, the corresponding 
and may be neglected regardless of its own particular 

magnitude. 

Since we are particularly interested in the significance of 
only [Am] and (Akl for the time being, we may examine only 
the real part, R[AZ] , of [A 21 initially. 


We may reason that a suitable criterion for neglecting Am^^ 
and Ak^j is that 

|{e}T[D](IAkl - a[Aml)[D](e}jl < e (4-26) 

where e ^ denotes the j column of the identity matrix, 

[l], and e<<l. That is, if we define 

[A^l] s [d](r[A2))[d], 

equation 4-26 may be stated 


Following similar logic to the derivation shown in Reference 
I 4 1 , we focus our attention on the resonant frequencies 
where we assume that 0 * < Wj? 5 The criterion 

for neglecting Am^^ and Ak^^ is thus found to be 



Vw J 
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To simplify the notation, we may define 
(Dj/(di • e > 1 (frequency ratio) 
l/2i;^ » (dynamic amplification) 

Then, for light damping (?j<<l) whenever (0^ - 1)^ >> 
equation 4-27 reduces to 

2 

As a practical matter, we might expect a value of e = 0.2 
to work. 

It is of interest to note from Reference [ 4 ] that the 
companion criterion for neglecting is 

I ^i 

Ac^j/w . ./ — 5 < e (4-29) 

1-3 ^ (g2 - 1) 

To conclude this section, we recall that our original objective 
was to band the matrices [m] and [k] for computational purposes. 
Since the rows and columns of these matrices are ordered by 
frequency, it is evident that the frequency ratio, B * 
which appears in the denominator of equation 4-28, acts to 
attenuate the importance of end Ak^^j as one gets 

further away from the diagonal. 



Sample Calculations 


The banding criterion derived in the previous seciton may 
be applied to some of the examples worked in Section 3. 

In particular, there were two runs made with banded [m] and 
(kl matrices. One converged and the other did not. 
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The first example corresponds to Example 2, Run 8, of 
Section 3.2. This run did converge in the sense that numerical 
computations stabilized after achieving a significantly im- 
proved fit of the data. The example consists of a two degree- 
of- freedom system with 10% damping in each moder where 


.0618 

-.0104 

-.0104 

1.0132 

‘ 1.0 

o.o“ 

, 0.0 

1.0. 


Since [®Ic] and [®m] are both diagonal (i.e. off-diagonal elements) 
are identically zero) in the prior model, it follows that 


Aki2 * -.0104 


'^">12 * ® 


It is also given in this example that 


- .0618 


coj » 1.0132 


6^ - 16.39 


corresponding to the diagonal terms of [h] , and 
0. 1 » 5.0 

We wish to determine whether in fact, equation 4-29 is satisfied. 


' ^1 
( 0 ^ - 1 ) 


.097 
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} 


which is less than the tentative value of e « 0.2 suggested 
earlier, i.e« the criterion worked. 

The second example corresponds to Example 3, Run 2, of 
Section 3.3. In this case, numerical computations did 
not converge when the (1,3) elements of Im] and Ck] were 
excluded from the set of parameters being estimated. The 
following data from Example 3 apply: 


Akj^3 » -3.35 
An»2^3 * -.136 
« 7.35 
» 118,5 
*0.05 


8^ - 16.12 
» 10.0 


We wish to determine again whether equation 4-28 is satisfied. 



» 0.481 


which is greater than the tentative value of e « 0.2 suggested 
on Section 4.2, i.e. the criterion worked again. 

The suggested value of e * 0.2 is only a first guess at 
establishing a working criterion for excluding off-diagonal 
elements of {m] and [k] from the vector of parauneters to be 
estimated in Phase I. In fact, this particular value was 
selected with the objective of achieving convergence, not 
necessarily with the objective of achieving a particularly 


i 

■1 


■J 


j5 


I 
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close fit of the data, or with the objective of achieving 
revised parameter estimates within any specified bounds. 
These latter two considerations may very well dictate a 
cutoff value for epsilon (e) of somewhat less than e * 0.2, 
It is reasonable to expect that the percent error Introduced 
in the estimates will be on the order of 100 x e, e.g* 20% 
for e = 0.2. As a matter of fact, in the first example of 
this section where an error value of .097 was calculated, 
response derivation of the revised analytical model from 
the "exact" solution was approximately 10%, 
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5. 


DEFINITXOH OF SUBMATRICES 


Phase II of the procedure described In Section 3 uses submatrix 
scaling coefficients as parameters of [M] and [K] to be esti- 
mated. No discussion was included with regard to how the sub- 
matrices are defined. Initially^ it was presumed that the anal- 
yst responsible for generating [H] and [K] would be able to 
define approximate submatrices, since the coordinate system 
associ^*‘ed with these matrices of the "dynauaic model’* does 
refer to physical displacements on the structure itself. This 
may still be one alternative. However, sufficient question 
was raised about the practical feasibility of doing this, that 
other alternatives must at least be considered. 

5.1 Sub structuring Basis 

It is well known that in linear structural analysis, mass 
ans stiffness matrices are generated by a process of super- 
position, When assembling these matrices from finite element 
contributions, the element mass and stiffness matrices are 
mapped into the global mass and stiffness matrices, in a 
similar manner, when substructuring techniques are employed, 
the substructure matrices are mapped into the complete struc- 
ture (global) system mass and stiffness matrices. 

This procedure suggests that each of these substructure 
matrices, after it has been mapped into global matrix elements, 
constitutes an appropriate submatrix for use in Phase II 
parameter estimation. In effect, if each substructure matrix 
is associated with a scaling parameter, the estimator will 
scale the overall stiffness and mass of the substructure up or 
down as requried to obtain a better fit of the data. Focus, 
with regard to how large a portion of the structure to look 
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at in each case, will of course depend on the relative size 
of the substructure compared to the system as a whole. In 
theory, it should be possible to subdivide the structural 
system however one chooses. In a practical sense, however, 
this option may be precluded by having to select the substruc- 
turing for other reasons, and perhaps by having not kept the 
intermediate results necessary to define the required sub- 
matrices. 

In conclusion, while the substructuring approach is theoretically 
sound and quite straightforward, it may sometimes be practically 
difficult to implement, particularly if one comes along after 
all of the modeling and condensation has been completed, and 
wishes to use intermediate results which are difficult to 
access, or not available at all. In this case, still another 

5.2 Orthogonal Modes Basis 

A third alternative has been postulated. Considering that all 
one may have to work with is IM] and [K] in their final form, 
the question is, is there any way these system matrices might 
be manipulated to derive appropriate submatrices for Phase II 
estimation. 

The use of orthogonal modes provides a basis for doing so. 

It is well known that the natural dynamic mode shapes of a 
structure can be used to generate a modal expansion of both 
the mass and stiffness matrices. The form of such an expan- 
sion is as follows: 

[HJ « Z [MJ {<}>} [M] (5-1) 

i i i 
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(5-2) 


[Kl - £ ojJ iVL] (Ml 

Assuming the orthogonality condition [M] * [ I ]. The 

implication here is that each submatriXf say 

[M]^ * (Ml [Ml, or 

(Kl^ * wj [M3^ 

could be associated with a scaling parameter, 3^ or a^, such 
that 


[Ml - £ [Ml^ 
[Kl * £ <»i (Kli 


In case a partial set of modes is used, we may irrite 


[Ml 


[M] 

+ E 5. 
i ^ 

[M]^ 

(5-3) 

[Kl 

m 

[Kl 

+ I a. 
1 * 


(5-4) 


where [M] and [Kl are defined in such a way as to complete 
[Ml and [Kl. 

This approach would conceivably work for the problem at hauui, 
with the a* 3 and 3’s initially set to unity. The estimator 
would adjust them according to the procedure described for 
Phase II. The only drawback to this particular approach is 
that it may not provide much insight into the meaning of 

adjustments made to elements of [M] and [Kl . 
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An alternative choice of modes for defining submatrices of 
£K] is one which satisfies the orthognality condition: 

{«}* [K1 {«), - 0 

(KI - K* 

Suppose that 

{6}^ - [K]"^ (5-6) 

where is a set of forces which are in equilibrium. That 

is, given {f}^ along with [K] , we can evaluate and K^. 

If we generate "modal deflections” 


<5-5a) 

(5-5b) 


{4»>i - i {&), 

i 


then 


[<(»] [Kl [<J>] * [ I ] 


and 




Thus, we may define 


[KJ - mM - i:[K], 

i ^ 


(5-7) 


(5-8) 


(5-9) 


(5-10) 


74 - 


where 


[K]^ * [KH<|>}^ [Kl 

KJ = (5-11) 

By definition, then 

IKl - (K) + r o^lK]^ (5-12) 

where the a's are the parameters associated with [K] to be 
estimated, and their initial values are 

» 1/K* (5-13) 

Our objective is to define force vectors such that (6}j*0 

i.e., the virtual work done by one set of forces, on the 
displacements caused by another set of forces, is zero. In- 
tuitively this will be the case, for example, when each set 
of forces is in equilibriiira, and applied to a different part 
of the structure. This approach should provide a means of 
defining submatrices associated, therefore, with different 
parts of a structure. 
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CONCLUSIONS AND RECOMMENDATIONS 


The following twol. subsections itemize specific conclusions and 
recommendations resulting from the study to date* 

6.1 Conclusions 

The formulation originally proposed for Phase I of the estima- 
tion procedure has been modified in the sense that both mass 
and stiffness parameters# as well as damping parameters, are 
now estimated. Originally# it was reasoned that by normalizing 
the modal vectors of the original "dynamic model" to yield 
an identity matrix for the mass matrix of the prior "analytical 
model#" that only elements of the (generalized) stiffness matrix 
of the "analytical" model would have to be estimated. This 
reasoning proved to be incorrect# and the necessary revisions 
to the procedure were made. 

Based on this modified formulation# a general computer code 
was developed. It is presently separated into two parts# one 
pertaining to the Phase l estimation procedure which refines 
a modal representation of the dynamic model# and the other 
pertaining to the Phase II estimation procedure which takes 
the refined modal representation and revises the mass and 
stiffness matrices of the original dynamic model. The two 
separate programs will be linked by a computer data file. 

Numerous example problems have been run to test the operation 
of computational algorithms. While some of these examples are 
not realistic in the sense that they do not necessarily repre- 
sent anticipated applications to real problems# they do serve 
to identify potential problem areas# and have been selected to 
test the methodology under extreme conditions. Thus, while 
the basic methodology has been successfully demonstrated, 
several related problem areas have also been illuminated. 
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Generally speaking, these problem areas can all be identified 
with one basic problem - that of determining a sufficiently 
good prior model frop which to begin the estimation procedure. 

A "good" prior model must first of all embody a " proper config- 
uration ". That is, it must contain the proper elements. If 
the modal representation of the prior "analytical model" is 
deficient in the sense that critical off diagonal elements 
(modal coupling elements) are ignored, then computations may 
not converge. If the critical elements are included, but in- 
itial estimates of their values are too far off, computations 
may still not converge. Both of these conditions tend to be 
aggravated by improper estimates of the uncertainty of the 
various parameters. For example, if the initial estimates 
of these parameters are greatly in error, while corresponding 
uncertainties are small by comparison, the two sets of input 
are injonsistent with each other, relative to the target model. 

As with any new tool, experience will be required to become 
proficient in its use. However, improvements to the method- 
ology are suggested by the experience gained so far. In 
particular, the analysis presented in Section 4 has the poten- 
tial of significantly improving (1) the configuration of the 
prior model, (2) the initial estimates of the prior model 
parameters, and (3) estimates of the uncertainties corresponding 
to these initial parameter estimates. 

The following subsection itemizes specific recommendations 
for subsequent investigation and implementation. 

6.2 Recoramendations 

Our first recommendation is to implement the first-order cor- 
rection procedure to the prior "analytic model described in 
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Section 4.1. This should take advantage of test data not 
otherwise used explicitly in the estimation procedure# to 
produce a significantly improved prior model. The efficacy of 
this approach should be easy to demonstrate by application to 
some of the analytic examples presented in Section 3. This 
approach should provide# in addition to improved estimates of 
initial parameter values# a quantitative basis for evaluating 
parameter uncertainties based in part on measured test data. 

Secondly# it is recommended that the banding criteria developed 
in Section 4.2 be investigated further as a means of helping 
to establish a proper model configuration. These criteria 
apply to the damping parameters, as well as the mass and stiff- 
ness parameters. 

In order to maximize the utility of the procedure, it is recom- 
mended that a nested sequential estimation capability be in- 
vestigated for implementation in Phase I. The intent here is 
to estimate parameters associated with a limited frequency band 
which moves incrementally across the frequency spectrum in the 
inner loop of the estimator# and then if necessary# sweep the 
frequency spectrum again with additional test data. This would 
in general expand the capacity for the number of parameters to 
be estimated. 

It is recommended that a method such as that described in 
Section 5.2 be implemented to define appropriate submatrices 
for use in Phase II. This would provide a straightforward 
means of localizing adjustments to the dynamic model# given 
only the mass and stiffness matrices of that model. 

Finally# it is recommended that additional analytical examples 
be investigated as a means of gaining experience with the 
modified programs prior to their application to real problems. 
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This experience should have a pay-o£f In helping to establish 
guidelines for all propsective users. 
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NOMENCLATURE 


The following symbols are used in this report; 


A 

a 

B 

b 

C 

c 

f(t) 

g(t) 

H(fl) 

H(£l) 

I 

M 

m 

P 

P 

q 

R 

r 

T 

t 

u 

U 

1C 

Y 

y 

z 

X 


» Real portion of H(n) • dynamic model coordinates 
a Real portion H(Q) - analytic model coordinates 
=s Imaginary portion H{£2) - dynamic model coordinate 
Imaginary portion H((2) - analytic model coordinate 
- Damping parameter - dynamic model 
» Generalized damping parameter - analytic model 
» Forcing function in dynamic model coordinates 
» Scalar function of time 

^ Transfer function - dynamic model coordinates 
» Transfer function - analytic model coordinates 
= Identity 

“ Mass parameter - dynamic model coordinates 
= Generalized mass parameter - analytic model 
= Force distribution parameter - dynamic model coordinates 
=s Force distribution parameter - analytic model coordinates 
» Coordinate of modal model, time domain 
= Effective state variable (equation 2-1) 

= State variable 

» Sensitivity of observation quantity to change in 
state variable 

a time 

a Coordinate of dynamic model, frequency domain 
a Coordinate of dynamic model, time domain 
a Coordinate of analytic model, time domain 
a observation (i.e., measured) quantity 
a Effective observation quantity (equation 2-2) 
a Variable of convenience (equations 2-10 and 2-27) 

a Coordinate of analytical model, freq^iency domain 
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I 


a 

6 

e 

n 

X 

C 

z 

a" 

n 

u 

A 

X 

[ 1 
{ } 


Scaling par 2 Lineters for dynamic model mass and 
stiffness matrices 

Kronecher delta function, also displacement vector 

Random error 

Damping factor 

Variable of convenience 

Eigenvalue (rad /sec ) 

Modal damping parameter - modal model 
Summation 

Variances or covariances 
Modal deflection - dynamic model 
Modal deflection - analytic model 
Excitation frequency (rad/sec) 

Natural frequency { rad/sec) 

Square root of -1 

Matrix 

Vector 


Supplementary Symbols 


AM 

AK 

AX 

A(J) 

An 

Am 

Ak 

D 

D(iJJ) 

r 

a 

0 


First order correction to (®MJ 

First order correction to [®K] 

Difference between prior model eigenvalues and 
measured frequencies squared 

Difference between prior model eigenvectors and 
measured mode shapes. 

Modal contribution (prior model modes) to A<^ 

First order correction to [®m] 

First order correction to [*kl 
Denominator 

Diagonal scaling matrix (complex) 

Scaled coordinate of analytical model# frequency domain 
Scaled coordinateof analytical model# time domain 
Frequency ratio; also used to denote scaling parameter 
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1 


I 


Subscript 


e 

g 

h 

i 

j 

k 

I 

m 

n 

o 

P 

q 

r 

3 


Desired value 

Summation index 

Index for state variables 

Index for dynamic model coordinates - response point 

Summation index, index for excitation frequencies, 
or mode shapes 

Index for dynamic model coordinates > driven point 
Index for "observation” vector 

Index for analytic model coordinates, summation index 
Summation index 
Summation index 
Prior model 

Index for analytic model coordinates, summation index 
Index of mass and stiffness matrices 
Index of mass and stiffness matrices 


Superscript (either right# left, or above) 


I 

N 

R 

T 

O 

it 


Imaginary component of complex quantity 
Total number of scaling parameters 
Real component of complex quantity 
Transpose of 

Original dynamic model, unvarying 
Complex conjugate 
Differentiation 

Non-vairying portion of matrix; also used 
normalized coordinates in frequency domain 

Analytic model 

Approximate Value (except for 


to denote 


Analytical model in scaled coordinate 
system 
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APPENDIX A 

DESCRIPTION OF COMPUTER PROGRAMS 

The parameter identification methodology and corresponding 
computer programs have been set up as a two step operation. 
First, a new set of generalized mass and stiffness matrices* 
are developed for the analytical model. Then these matrices are 
used as "observation” data to develop new scaling parameters 
for the dynamic model mass and stiffness matrices. 

The two computer programs are not automatically coupled. Cur- 
rently, the user must transfer the results from Phase I to 
cards to be used as input to Phase II. Eventually the identi- 
fication programs will be configured so that disk files can be 
used as the transfer medium removing much of the work load from 
the user. Each program would still function independently, how- 
ever. The advantages to keeping each phase independent are 

e It may take several runs before Phase I produces 
a set of inodes and frequencies acceptable to the 
engineer from his knowledge of the structure. 

e The user should evaluate the Phase I results 
before proceeding to Phase II to insure the new 
model is meaningful. 

e The computer programs are simpler and the core usage 
more efficient. 

Both programs were originally dimensioned for relatively small 
problems to minimize computer costs during the early check-out 


The generalized mass and stiffness matrices for the dynamic 
model are the mass and stiffness matrices for the analytic 
model. 
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4 I 


and investigation phase. They will be redimensioned in the 
final versions of the programs to take advantage of available 
computer memory core when running the practical demonstration 
problems. 

A-i Program ESTIMA (Phase I) 

The Phase I procedures have been subdivided in two independent 
operations which are executed alternately and repeated until 
convergence to the test response is obtained. Phase I*-A iter- 
ates on the mode shpaes and natural frequencies until no further 
improvement can be obtained. Phase I-B then iterates on the 
modal damping to further improve the fit. When no further 
improvement is obtained here. Phase I-A is repeated. These 
two operations alternate until the optimum combination of modal 
parameters (frequencies, mode shapes, damping) is obtained. 

The equations needed to perform these operations have been 
programmed for a CDC 6600. These have been broken down into 
a main program plus nine subroutines. A flow chart of the 
main program is shown in Figure 1. As of this date the portion 
of the computer program specifically related to Phase I-B oper- 
ations has not been entirely verified. The example problems 
perfozrmed to date have all been designed to develop a greater 
understanding of the Phase I-A methodology. They have been 
set up as if the damping ratios were exactly known prior to 
the identification process. 

DISRES Computes transfer functions for analytic 

model coordinates and dynamic model degrees 
of freedom. Computes frequency response for 
K same two coordinate systems. 
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FIGURE 1 . FLOW CHART FOR PROGRAM ESTIMA 



Q START J 


READ ALL INPUT DATA RELATING TO PROGRAM 
CONTROL AND PRIOR MODEL 




/ 


WRITE ALL INPUT DATA RELATING TO PROGRAM 
CONTROL AND PRIOR MODEL ON PRINTER 



r 

^ “ 1 

1 / READ TEST DATA FOR THE INITIAL oHMATION / 1 

1 Z ^ 1 

1 

1 


1 /mm test data for the initial estimation/ i 

I y OH PRINTER . / • 

1 

J 


340 


SET ICHK TO TRUE OR FALSE, ICHK CONTROLS 
CERTAIN OPERATIONS THAT ARE ONLY PERFORMED 
ON FIRST ITERATION 


SET UP THE INITIAL M MATRIX EQUAL TO 
IDENTIFY. M IS THE MODAL MATRIX OF THE 

analytic system [nm by NM] 
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?T-I J|T 





FLOW CHART FOR PROGRAM ESTIMA (Continued) 



HOT EXCEEDED.^tf«?T rvcLF^Vw EXCEEDED 
i S^LIMIT NCLLMT^ 


Is THE ycc 
’RINT CONTROL, JPRINT>^ 
^ TRUE ? 


WRITE INTERMEDIATE 
MATRICES ON THE PRINTER 



FORM THE EFFECTIVE "OBSERVATION" VECTOR 
Y{J) - Uy^syCJ) - 


WRITE THE OBSERVATION VECTOR ON 
THE PRINTER 


T7-I 2IT 




aON CHART FOR PROGRAN ESTIMA (Continued) 



HAS RMS ^ 
ERROR INCREASED FROM 
PREVIOUS ITERATION ? 


RESET SIZE CONTROLLER TO 
1/4 OF PREVIOUS STEP 


TRANSFER NEW ANALYTICAL RESPONSE TO U FOR 
STORAGE 


' SUBROUTINE SENSK OR SENSO 


FORM THE SENSITIVITY MATRIX [T] 


^ IS THIS ^ 
THE FIRST CYCLE 


TRANSFER THE NEW PARA- 
METERS, R, TO PRIOR 
PARAMETERS, RP 


Id 


FORM THE VECTOR, RP, OF THE PRIOR 
PARAMETERS 
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FLOW CHART FOR PROGRAM ESTIMA (Continued) 
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FLOW CHART FOR PROGRAM ESTIMA (Continued) 



TT-tJ*T 






C>''y 


FLOW CHART FOR PROGRAM ESTIMA (Continued) 
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SENSK 


SENSD 

MOUSE 

GIVHO 

INVERT 

INVECC 

MATMUL 

MATOUT 


Computes the sensitivity matrix for analytic 
model mass and stiffness matrices (i.e., 
the sensitivity of the frequency response with 
respect to variations in the analytical model 
mass and stiffness matrix) . 

Computes the sensitivity matrix for modal 
damping matrix. 

Estimates a new analytic model stiffness matrix. 
Also develops a new covariance matrix when each 
set of iteration cycles is complete. Based on 
same equations as computer program MOUSE [2] • 

Uses Householders method to reduce a real 
symmetric matrix to tridiagonal form. Iso- 
lates eigenvalues and eigenvectors* Adapted 
subroutine from program MOUSE [2]. 

Inverts a real symmetric matrix. Adapted 
subroutine from program MOUSE [2] . 

Inverts a complex symmetric matrix. 

Multiplies tow conformable matrices. Adapted 
subroutine from program MOUSE [2]. 

Prints all non-zero elements of a matrix. 

Adapted subroutine from program MOUSE [2] . 
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INPUT DATA FOR PROGRAM ESTIMA 


Card One 

Cards 1 to 3 provide information on input<*output requirements^ 
on the convergence criteria/ and the elements to be estimated. 

IPRINT, JPRINT 

MMW OTM* wna 4mrnmm smmpv « 

LIO / LlO FORMAT 

where 

XPRXNT is a flag which controls printing of the 

modal orthogonality check. 

s»T print the entire orthogonality check 

5*P print only bad elements of orthogon- 
ality check 

JPRINT is a flag which controls printing of the 

various intermediate matrices 

=*T print intermediate matrices 

ap do not print intermediate matrices 

Card Two 


NCLLMT/ 

MB 


15 

15 

FORMAT 


where 

NCLLMT is the maximum number of iterations 

allowed 

NB is the bandwidth within which new ele- 

ments are to be estimated. 
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Card Three 


CONLMTf CONLM2 


FIO.O , PIO.O FORMAT 


where 

CONLMT is a convergence criterion (ratio between 

calculated and observed response that 
will constitute success) . 

COMLMT is a second convergence criterion (ratio 

between two successive calculations 
that will constitute success) . 


Card Four 

Cards 4 to 11 provide data on the prior model. 


NM , NC 


15 ,15 FORMAT 


where 

NM is the number of modes being used. 

NC is the number of degrees of freedom in 

the u coordinate system. 

Card Five 


W(I) , I » 1 NM 


8F10.0 


FORMAT 








where 


NDMPFL la the control flag which specifies the 

type of damping information to be read* 

sQ read the critical damping ratios for 
the NM modes of the prior model 

*1 read the full NM by NH damping matrix 
for the prior model 


Card Seven (if NDMPFL » 0) 


ZET(I), I • 1, NM 


8F10.0 

FORMAT 


where 


ZET(I) are the critical damping ratios for the 

NM modes of the prior model being used 
for this estimation. 


Card Seven (if NDMPFL * 1) 


(ETA(I,J), I - 1, NM) , J « 1, NM 


8F10.0 FORMAT 


where 

ETA (If J) are the NM times NM elements of the modal 
damping matrix for the prior model being 
used for this estimation. 


Card Eight 


(PHI(I,J), I • 1, NO, J ■ 1, NM 


8P10.0 FORMAT 
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where 




are the elements of the modal matrix 
for the NM modes being used for this 
estimation. 


Card Nine 


SRPRP(I) , I » 1, NP 


SFIO.O 

FORMAT 


where 


SRPRF(I) are the initial variances of the elements 
of the generalized mass matrix that are 
to be estimated. The program determines 
NP based on the band width specified 
and the number of modes. Only the diagonal 
and upper right elements are estimated. 


Card Ten 


SRPRP(I), I e NP+1, 2*NP 


SFIO.O 

FORMAT 


where 


SRPRP(l) are the initial variances of the elements 
of the generalized stiffness matrix 
that are to be estimated. These elements 
correspond to the elements of the mass 
matrix being estimated. 


Card Eleven 


SRFRPD(I) , I - 1, NP 


SFIO.O 

FORMAT 










where 


SRPRPD(I) are the initial variances of the generalized 
damping matrix* These elements corres- 
pond to the elements of the mass matrix 
being estimated. 


Card Twelve 


Cards 12 to 14 provide the observation data (l.e,, test data). 
Successive sets of observation data may be processed. 


NF 


15 

FORMAT 


where 

NF is the number of excitation frequencies 

for which observation data is being read* 


Card Thirteen 


FQ(I) , 1 - 1, NP 
8F10.0 

FORMAT 

where 

PQ(I) 

are the observation frequencies (Hz) . 


Card Fourteen (a block of cards) 

This block of cards must be provided for each of the NF obser- 
vation frequencies. L is the index of the observation fre- 
quencies. 
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1 





card I4-A 


NCT(L), NCP(L) 


15 « 15 


FORMAT 


is the number of coordinates with obser*^ 
vation data. 

is the number of coordinates being forced 
(i.e.f with shakers). 


Card 14-B 


1CT(I,L) ,1-1, NCT(L) 
1615 


where 


FORMAT 


ICT(I,L) are the locations (i.e., the degrees 
of freedom) of the observation data 
(i.e., the measured response). 


Card 14-C 


ICP(I.L) ,1-1, NCP(L) 


1615 

FORMAT 


where 


1CT(1,L) are the locations (i.e., the degrees of 
freedom) of the coordinates being forced 








Card 14*0 


PR(l.L) Z - 1, NCF(L) 


8F10.0 


FORMAT 


where 


PR(I,L) 

are the real components of the excitation 
forces. 

Card 14*E 



P1(I,L) ,1-1, NCP(L) 

8P10. 


FORMAT 


PI{I,Ii) are the imaginary components of the 
excitation forces. 


Card 14-F 


UTEST(I), I - 1, MO 

- 

8F10.0 

FORMAT 


where 


UTEST(Z) are the observed responses arranged as 
followed: 

UTEST(l) » location 1, freq, 1 
UTEST(2) - location 1, freq. 2 

through all observation frequencies. 
Repeat for all locations. 







Card 14-G 


SEE (I) , 1 - 1, NO 


8F10. 


SEE (I) 


are the variances of the observed 
response arranged as described jibove. 

# . 
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A-2 Program ESTIMB (Phase II) 


' W 


ESTIMB, the Phase II, minimum variance estimator consists of 
a main program plus five subroutines. A flow chart of the 
main program is shown in Figure 2. The five subroutines ares 

ESTIMB Main program. Estimates a new dynamic model 

(mass and stiffness matrices) from the refined 
parameters of the analytical model (from Phase I). 

MOUSE Estimates new parameters (in this case the 

Oj) . Provides the covariance matrix of the 
new parameters. Identical to the MOUSE sub- 
routine used for Phase I. 

MMULRR Performs matrix multiplication of two real 
matrices without destroying either one. 

MATOUT Prints all non-zero terms of a matrix with 

pages, page headings, and matrix identifica- 
tion. Identical to the MATOUT subroutine 
used for Phase I. 

MAT0U2 Prints all non- zero terms of a matrix without 

pages, page headings, or a matrix identification. 

INVERT Inverts a real symmetric matrix. Identical 
to the INVERT subroutine used for Phase I. 

Uses the Choleski SDS decomposition method. 

Input Data for Prograun ESTIMB 

Card One 

Cards 1 to S provide information on input-output requirements, 
on the convergency criteria, and on the number of elements 
to be estimated. 
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FIGURE 2. FLOW CHART FOR PROGRAM ESTIMB 


Z 


START ^ 


READ ALL INPUT DATA RELATING TO 
PROGRAM CONTROL 


I 


7 


SET A NUMBER OF CONTROL PARAMETERS 




READ THE DYNAMIC MODEL MASS AND 
STIFFNESS MATRICES, [M] 4 [K], 1 
MATRICES FOR ERRONEOUS DATA 




^ SUBROUTINE MATOU2 

/ WRITE THE DYNAMIC MODEL MASS AND / 

/ STIFFJ^ESS MATRICES ON THE PRINTER / 



PUCE MASS MATRICES [M]^ ON TAPE 21 , 
STIFFNESS MATRICES [K]^ ON TAPE 22 


Z 




READ AND WRITE ON PRINTER COVARIANCES 
OF PRIOR SCALING PARAME F 5 



7 


TT-llIT 
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FLOU CHART FOR PROGRAM ESTIHB (Continued) 


L 


1 


SUBROUTINE MAT0U2 


READ AND WRITE ON PRINTER PRIOR MODE 
SHAPES 


I 


7 


SUBROUTINE MAT0U2 


/read AND WRITE ON PRINTER OBSERVA- 7 
/ TIONS - I.e., ELEMENTS OF GENERALIZED/ 
/ MASS AND STIFFNESS MATRICES; / 

/ [mj + M from phase I / 

I 

▼ SUBROUTINE MATOUT 

Z READ AND WRITE UPPER TRIANGUUR / 
HALF OF COVARIANCE MATRIX OF / 
OBSERVATIONS FROM PHASE I / 


I 


WRITE COVARIANCE MATRIX ON TAPE 
FOR INTERNAL STORAGE 




FLOW CHART FOR PROGRAM ESTIMB (Continued) 
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FLOW CHART FOR PROGRAM ESTIM8 (Continued) 
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FLOW CHART FOR PROGRAM ESTIMB (Continued) 
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JPRINT 



LIO 


FORMAT 

where 

JPRINT 

is a flag which controls printing of 
intermediate operations , 



^ print intermediate matrices. 

do not print intermediate matrices. 

Card Two 


■* 

NCYLLT 



15 


FORMAT 

where 

NCYLLT 

is the maximum number of iterations 
allowed. 

Card Three 


CONLMT, 

C0NLH2 


PIO.O , 

FIO.O 

FORMAT 


where 


CONLMT is a convergence criterion (ratio between 

calculated and observed values that will 
constitute success) for stiffness matrix, 

C0NLM2 is a second convergence criterion (ratio 

between two successive calculated values 
that will constitute success) for stiff- 
ness matrices. 
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I 


Card Four 


CONLM3, 

C0NLM4 


FIO.O , 

FIO.O 

FORMAT 

where 

C0NLM3 

is a convergence criterion (ratio between 
calculated and observed values that will 
constitute success) for mass matrix. . 


C0NLM4 

is a second convergence criterion (ratio 
between two successive calcualted values 
that will constitute success) for mass 
matrix. 

Card Five 



NK, NM, 

NS 


15, 15, 

15 

FORMAT 

where 

NK 

is the number of K-matrix scaling para- 
meters to be estimated. 


NN 

is the number of M-matrix scaling para- 
meters to be estimated. 


NS 

is the size (degrees-of-freedom) of the 
dynamic model. 


Card Six 


Cards 6 and 7 are repeated NK times. Cards 6 to 12 provide 
information on the prior model. 
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N 

15 

FORMAT 

where 


N 

is the index number of the K-matrix 


portion to be read next. 

Card Seven 


JJ, KK, K(JJ,KK) 


4(I5,XS,F10.0) 

FORMAT 


where 

JJ 
KK 

K(JJ,RK) 

Values of JJ, KK, and K(JJ,KK) are read until a value of 0 
Is read for JJ. Use as many cards as necessary with 4 elements 
per card. 

Repeat cards € and 7 until all portions of the prior K-matrix 
are read in. The program will continue to read K-matrix blocks 
until a value of 0 is read for M. Only the upper right-hand 
elements of the K-matrices must be read. 

Card Eight 

Cards 8 and 9 are repeated NH times. 


is the first index of the element. 

is the second index of the element. 
KX must 3^ JJ. 

is the JJ,KX element of the portion 
of K-matrix being read. 
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15 


FORMAT 

where 



N 

Card Nine 

is the index number of the 
portion to be read next. 

M-matrix 

JJ, KK, M(JJ,RK) 

4(15, 15, FlO.O) 


FORMAT 


where 

JJ is the first index of the element. 

KK is the second index of the element. 

KK must ^ JJ. 

H(JJrKK) is the JJ,KK element of the portion of 
the mass matrix being read. 

Values of JJ, KK, and H(JJ,KK) are read until a value of 0 is 
read for JJ. Use as many cards as necessary with 4 elements 
per card. 

Repeat cards 8 euui 9 iintll all portions of the prior H-matrix 
are read in. The program will continue to read M-matirx blocks 
until a value of 0 is read for 19. Only the upper right-hand 
elements of the H-ma trices must be read* 
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Card Ten 


SRPRP(X) , X « HP 


8FX0. 


FORMAT 


SRPRP 


are the variances of the original scaling 
parameters. HP • NM + NK - 2. 


Use as many cards as necessary to read all of the variances. 

A maximum of 20 values are presently allowed^ 10 for K<*matrlx 
p^arameterSf 10 for H-matrix parameters. 

Card Eleven 


ND 


15 

FORMAT 


where 



is the number of modes of the dynamic 
model that are being used. 

Card Twelve 



PH1P(J,I) , I - 1, HS 


8F10. 


FORMAT 


where 


PHlP(Jfl) is the element of the mode 













Use as many cards as necessary to complete each mode shape. 
Repeat Card 12 for each mode shape. A maximum of 30 modes with 
30 elements per mode may be read. 


Card Thirteen 


Cards 13 to 15 provide the "test” data. 



where 

10 is the 1 subscript. 

JO is the j subscript. 

KO is the if j element observation general- 

ized stiffness matrix [k] • 

HO is the if j element of the observation 

generalized mass matrix r [ml . 

Repeat Card 14 as often as necessary to read all of the "test* 
data. A maximum of 30 elements may be used. 
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Card Fifteen 


JJ, RK, S£E(JJ,RR) 


4(15,15, E10.4) 

FOPMAT 


where 

JJ is the first index of the elemer.*-. 

KK is the second index of the element. 

KR must ^ JJ« 

SEE(JJfKK) is the JJ«KK element of the covariance 
matrix of the "observation** data. 

Only read the upper right-hand elements of the covariance 
matrix. 
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